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Abstract 

In this article we present new results for the pricing of arithmetic Asian 
options within a Black-Scholes context. To derive these results we make 
extensive use of the local scale invariance that exists in the theory of con- 
tingent claim pricing. This allows us to derive, in a natural way, a simple 
PDE for the price of arithmetic Asians options. In the case of European 
average strike options, a proper choice of numeraire reduces the dimension 
of this PDE to one, leading to a PDE similar to the one derived by Rogers 
and Shi. We solve this PDE, finding a Laplace-transform representation 
for the price of average strike options, both seasoned and unseasoned. This 
extends the results of Geman and Yor, who discussed the case of average 
price options. Next we use symmetry arguments to show that prices of 
average strike and average price options can be expressed in terms of each 
other. Finally we show, again using symmetries, that plain vanilla options 
on stocks paying known cash dividends are closely related to arithmetic 
Asians, so that all the new techniques can be directly applied to this case. 
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1 Introduction. 



Options with a payoff- function which depends on the average of some underly- 
ing, a.k.a. Asian-type option have a multitude of applications in finance. They 
find applications for example in currency-based contracts, interest rates and 
commodities. In the following we will consider a setting where stock prices 
are modeled by geometric Brownian motions. Depending on the type of av- 
eraging the analytic price of such a contract is easy or diffic ult to c ompute. 
Geometric averaging leads to simple expressions for the prices |HN99b . Arith- 



metic averaging however is a highly non-trivial exercise and one has to rely 
on either approximations |TW91 |, pa rtially analytic [GY92] or numerical so- 
lutions | Cur94 , RS95 |. In Rcf. RS95 ] a simple PDE was derived and bounds 
on prices of an average price option were derived. These bounds have been 



improved in Ref. [ Lho 



In this article we provide an alternative approach to derive (partially) analytic 
solutions of Asian-type contracts with arithmetic averaging. Using the funda- 
mental notion of a local scale invariance [HN99a, HN99b| we derive a general 
PDE for a (European) Asian-type contract with arithmetic averaging. This re- 
sult is then linked to the PDE derived by Rogers and Shi [ RS95 1 and we proceed 
by solving the Laplace-transform of the solution of this PDE for the case of an 
average strike put, both for the unseasoned and seasoned case. This extends 
the results of Geman and Yor, which gave the solution for average price options. 
Next we show that the local scale invariance allows one to identify the average 
strike call and average price put by substitution of the proper parameters. This 
is a new result. Finally we consider the problem of a vanilla option on a stock 
which pays known cash dividends. Again the local scale symmetry allows one 
to relate the value of such a contract to that of arithmetic average options. 



2 Homogeneity and contingent claim pricing 



In previous papers |HN99a| , [HN99b1 we have shown that a fundamental property 



of any properly defined market of tradables^ is that the price of any claim de- 
pending on other tradables in the market should be a homogeneous^ function of 
degree one of these same tradables. This property is nothing but a consequence 
of the simple fact that prices of tradables are only defined with respect to each 



other. Let us review some of the content of Rcf. | HN99a |. Assume that we have 
a market of n + 1 basic tradables with prices (fj, = 0, . . . , n) at time t. The 
price of any tradable in this market with a payoff depending on the prices of 

Tradables are objects which are trivially self-financing: it doesn't cost nor yield money 
to keep a fixed amount of them. Examples are stocks and bonds. Note that money is not a 
tradable, unless the interest rate is zero. 

2 A function f(xa, ... ,x n ) is called homogeneous of degree r if f(axQ,... ,ax„) = 
a r f(xQ,... ,x n ). Homogeneous functions of degree r satisfy the following property (Euler): 

c f ax 
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these basic tradables should satisfy the following scaling symmetry: 



V(Xx,t) = XV(x,t) 
which automatically implies^] (Euler) 

V(x,t) = x^d Xfi V(x,t) 

where d/dx^ = d x . This is a universal property, independent of the choice of 
dynamics. We use this fundamental property to derive a general PDE, giving 
the price of such a claim in a world where the dynamics of the tradables are 
driven by k independent standard Brownian motions, as follows^] 

dx^{t) = (<r M (a;, t) ■ dW(t) + a^(x, t)di) , (no sum) 

Consistency requires that both <r M and are homogeneous functions of degree 
zero in the tradables, i.e. they should only depend on ratios of prices of trad- 
ables. Note that we do not specify the numeraire in terms of which the drift 
and volatility are expressed. This choice is irrelevant for the pricing problem, 
as we will see. Applying ltd to V(x,t) we get 

dV(x, t) = d Xfi V(x, t)dx^ + CV(x, t)dt 

where 

CV(x,t) = (d t + i(Tp(a;,i) • a y (x,t)x ll x u d Xil d x ^jV 

So, if V(x,t) solves CV — with the payoff at maturity as boundary condition 
V(x,T) — f(x), then we immediately have a replicating self-financing trading 
strategy because of the homogeneity property. We will drop the distinction 
between such derived and basic quantities and always refer to them as tradables. 
Note that we do not have to use any change of measure to arrive at this result, 
by keeping the symmetry explicit. Drifts are irrelevant for the derivation of 
the claim price. Only the requirement of uniqueness of the solution, i.e. no 
arbitrage, leads to constraints on the drifts terms if deterministic relations exist 



between the various tradables [HN99a 



2.1 Symmetries of the PDE 

The scale invariance of the claim price is inherited by the PDE via an invariance 
of the solutions of the PDE under a simultaneous shift of all volatility-functions 
by an arbitrary function X(x,t) 

a-fj,{x,t) — > (Tn(x,t) - X(x,t) (1) 

3 We make use of Einsteins summation convention: repeated indices in products are im- 
plicitly summed over, unless stated otherwise. 

4 Both the <r M and dW are vectors, the dot denotes an inner-product w.r.t. the k driving 
diffusions. 
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Indeed, if V is solves LV = 0, then it also solves 

(®t + ^(^/t(»>*) - A(ar,t)) • (<J„(a;,i) - A^^a^a^d^a^J V = 

This can easily be checked by noting that for homogeneous functions of degree 1 
we have 

x^d x ^d Xu V = 

This ensures that terms involving the A drop out of the PDE. (Note that this 
equation gives interesting relations between the various T's of the claim). From 
this it follows that V itself must be invariant under the substitution defined by 
Eq. [|. This corresponds to the freedom of choice of a numeraire. It just states 
that volatility is a relative concept. Price functions should not depend on the 
choice of a numeraire. 

2.2 The algorithm 

To price contingent claims we start out with a basic set of tradables. Using 
these tradables we may construct new, derived, tradables, whose price-process 
V depends upon the basic tradables. Of course, these new tradables should be 
solutions to the basic PDE, CV — 0. Their payoff functions serve as boundary 
conditions. (Note that prices of basic tradables trivially satisfy the PDE, by 
construction). If the derived tradables are constructed in this way, we can use 
them just like any other tradable. In particular, we can use them as underlying 
tradables, in terms of which the price of yet other derivative claims can be 
expressed (and so on...) In fact, this is a fundamental property that any correctly 
defined market should posses. It amounts to a proper choice of coordinates to 
describe the economy. 

The general approach to the pricing of a path-dependent claim in our formalism 
can be described as follows. 

1. The payoff is written in terms of tradable objects. 

2. A PDE is derived for the claim price with respect to these tradables. 

3. The PDE is solved. 

4. Possible consistency check: the solution should be invariant under the 
substitution Eq. [l] (numeraire independence). 

2.3 Generalized put-call symmetries 

As an example of the strength of this symmetry, and to show the natural em- 
bedding in our formalism, consider an economy with two tradables with prices 
denoted by x\^ and dynamics given by (i = 1, 2) 

dxi(t) = Xi{t)oi(x\,X2,t) ■ dW(t) + . . . (no sum) 
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It is easy to see that under certain conditions there should be a generalized 
put-call symmetry. Any claim with payoff f{xi,x^) at maturity and price 
V(x±,X2,t) should satisfy 

d t + ^\a(x 1 ,x 2 ,t)\ 2 x 2 1 d 2 x ^jV = 

where a{x\,X2,t) = ai{x\,X2,t) — a2(xi,X2,t). Homogeneity implies that it 
also solves 

d t + ^\a(x 1 ,x 2 ,t)\ 2 x 2 d 2 ^jV = Q 

Therefore, if \a(xi,X2,t)\ 2 = \<r(x 2 ,xi,t)\ 2 , this PDE can be rewritten as 

d t + ^\a( X 2,x 1 ,t)\ 2 x 2 d 2 ^V = 

and we see that V(x 2 , x-y, t) with payoff f(x 2 , x±) is a solution, too. This is noth- 
ing but a generalized put-call symmetry. In the first case X2 acts as numeraire, 
in the second case x\ takes over this role. The usual put-call symmetry follows 
if we take a constant a and let x\, x 2 represent a stock and a bond respectively. 

2.4 Lognormal asset prices 

In an economy with lognormal distributed asset-prices 

dx^(t) = x^a^t) ■ dW(t) + . . . (no sum) 

it is possible to write down a very elegant formula for European-type claims, as 
was shown in Ref. HN99a | 

V{x , ...,x n ,t) = J V(x (l)(z - Oo), x n <j>(z - 9 n ), T)d m z (2) 

with 



(3) = (^r° p HI 



z? 



The 9^ are m-dimensional vectors, which follow from a singular value decom- 
position of the covariance matrix £p„ of rank m < k: 
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3 Arithmetic Asians 



In this section we will consider the pricing of Arithmetic Asian options. Since 
this is the only type of Asian options that we will look at, we will omit the word 
'Arithmetic' in the sequel. Note that parts of this material already appeared in 



Ref. HN99b |. A fundamental building block in the construction of a European 
Asian option, expiring at time T, is a tradable which at time T represents the 
value of a stock at an earlier time s. But to define this, we must first agree 
how to translate value through time. For this, we need a reference asset. A 
convenient choice is to take a bond P(t,T) (or P(t) for short), which matures 
at time T, as reference. Then we can define 

Ys{t)= { m p ® 

In words, this is a portfolio where one starts out with a stock, and converts 
it into a bond at time t — s. It is trivially self-financing. To set the stage, 
we will assume that the interest rate has constant value r, as is usual in the 
Black-Scholes context (stochastic interest rates are much harder to handle, see 



Rcf. [HN99bQ. In that case the bond with maturity T has value e~ rtyT ~ 1 ^ when 
expressed in the currency in which it is nominated (say dollars) . Consequently, 
an amount e -r( - T_s ) of the tradable Y s will have a dollar value at time T which 
is equal to the dollar value of the stock at time s. We will also assume that 
the contracts are initiated at time t — 0, unless stated otherwise. Note that for 
t e [0, T] we have 

S(t)=Y T (t) 

P(Q) (3) 

3.1 Discretely sampled Asians 

The main reason for introducing the objects Y s (t) is that they constitute a 
natural basis of tradables in which prices and payoffs of arithmetic Asians can 
be expressed. For example, the payoff of a discretely sampled average price call 
(APC) can be written as 



E 



w(U)Y ti (T)-KP(T^J = (fc w (U)^-Kjp(T) 



where {ti} is a set of sample times, w(U) are corresponding weights, and K 
is the strike. We use the notation (•)+ for max(-,0). Observe that the payoff 
explicitly contains P(T) — 1$ to make it homogeneous of degree one in the 
tradables. In a similar way, the payoff of an average strike put (ASP) becomes 



6 



where A; is a (generalized) strike. In view of Eqs. ||, we see that both options 
are in fact instances of a more general discrete Asian option, which is defined 
by the payoff 



V({t i },w,T)= (^2 W ( ti )Y ti (T^j 



(4) 



3.2 Valuation by multiple integrals 



In this section we want to calculate the value of the generalized option defined 
by Eq. [|, at the time the contract is initiated, i.e. t = 0. This value is known as 
its unseasoned value. Without loss of generality, we will assume that there are 
N + 1 sample times, satisfying = to < ti < • ■ • < £iy— l < tjv = T. If we take 
the bond as numeraire, and assume that the stock price follows a lognormal 
price process 

dS(t) = aS(t)dW{t) + ■■■ 
then it is easy to derive that the tradables Y s (t) satisfy 

dY s (t) = l t<s aY s (t)dW(t) + --- 

where 1(< S is the indicator function. It equals one when t < s and zero other- 
wise. So the Y s (t) also follow lognormal price processes with a time dependent 



volatility and we can directly use the results of section 2.4. The first step is to 
calculate the variance-covariance matrix of the tradables 



/ l t<ti l t<tj dt = o- 2 min(ti,tj) 
Jo 



The rank of this matrix is N. A singular value decomposition = 9i ■ 9j is 
given by 

0o = ((V^O), 9i = (A,... , A, 0, .... 0) 

N N-i 

and the Pi are defined by 

Pi = o-y/U - U-i 

The value of the unseasoned option can now be written as an iV-dimensional 
integral (here z = (zi, . .. , z/v)) 



V({t l },w,0) = S(0) 



j d N z(^w{U)cj ) {z-e i )^ (5) 



where we used the fact that Y S (Q) = 5(0) for all < s < T. This is in fact a 
Feynman-Kac formula. 
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3.3 An interesting duality 

In this section, we will again calculate the unseasoned value of the option in 
Eq. |[ but this time we take the stock as numeraire. This corresponds to a shift 
of er in all volatility functions. We find that Y s (t) now satisfies 

dY s (t) = l t>s aY s (t)dW(t) + --- 

In this case the variance-covariance matrix becomes 

= cr 2 / UyuUytjdt — er 2 min(T - U,T - tj) 
Jo 

and a singular value decomposition Sjj = Bi ■ 9j is given by 

k = (0,... ,0,/Ji+i,... ,p N ), 6 N = (0,... ,0) 

i N 

So an alternative expression for the price of the option is 

V{{tih w > °) = S (°) J dNz (j2 W ^)<P( Z - Oi^J 

However, if we compare this result with Eq. 0, we see that it could also be 
interpreted as the value of an option with payoff 

In other words, two options which are related by the following substitution in 
their payoff 



Y t (T)~Y T -t(T) (6) 



have the same value at t = 0. We will call this T-duality (T from Time-reversal). 
It is a very interesting symmetry operation because, in view of Eqs. ^, we can 
use it to relate the values of unseasoned average strike and average price options. 
We will come back to this point in section 3.6. Note that Eq. ^ takes a simple 



form by virtue of the fact that we are working in a basis of tradables. 
3.4 A PDE approach 

In this section we consider a PDE approach to the pricing of Asian options. We 
derive a very general PDE, which can be related to the one that is usually found 
in the literature. Our PDE, however, has the advantage of being manifestly 
numeraire independent by virtue of the fact that it is expressed in a basis of 
tradables. It can be used to price both American and European style options, 



8 



but we will focus on the European case here. We will come back to American 
Asians in future work. The basic idea in the derivation of the PDE is, instead 
of introducing a tradable for each sample date, to introduce one new tradable 
S(t), which is a weighted sum over Y s (t) (obviously, a sum of tradables is again 
a tradable). This allows us to consider continuously sampled Asians in a proper 
way. Also 

S(t) = [ w(s)Y s (t)ds = <f>(t)S(t) + A{t)P(t) 
Jo 

where A(t) is proportional to the running average 

Ait> =l! w{s) w) d ' 

and 

4>{t) = w(s)ds < — > w(t) = -d t <t>{t), 4>{T) = 

Of course this approach also incorporates discretely sampled Asians. In that 
case, w(t) will be a sum of Dirac delta-functions and </>(t) will be a piecewise 
constant function, making jumps at sample dates. If we choose the bond P(t) 
as numeraire, and assume that S(t) satisfies 

dS(t) = aS{t)dW(t) + ■■■ 

then it is obvious that the new tradable S(t) satisfies 

dS(t) = cj>(t)aS(t)dW(t) + ■■■ 

This straightforwardly leads to the following PDE for options which depend on 
S, S and P 

(d t + \a 2 S 2 {ds + 4d § ) 2 ^V = (7) 

If we perform a change of variables in the PDE, eliminating S in favor of A, we 
find 

(dt + wpA + ^&^V = 

which is closer to the usual formulation. However, since A(t) does not corre- 
spond to the value of a tradable object, this form of the PDE looses the manifest 
symmetry. 

3.5 Analytical solutions 

In this section we derive a Laplace transform representation for prices of aver- 
age strike options in the case that sampling is continuous with an exponential 
weight function. This is in fact the counterpart of the calculation of Geman and 
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Yor [GY93] for the price of average price options, although they used entirely 
different methods to derive it. The results for average strike options turn out to 
be somewhat more involved, mainly because there is no simple relation between 
prices of seasoned and unseasoned options, while a simple relation does exist in 
the case of average price options, as we will see. Now, by definition, the payoff 
of an average strike option is defined in terms of S and S only, P does not 
appear. For example, the payoff of an ASP is given by 

(S(T) - kS(T))+ (8) 

Because of this fact, it is natural to choose S as numeraire. Dropping the 
derivative w.r.t. P, the PDE then reduces to 

d t + ^<T 2 (s-<f>s) 2 d 2 ^jv = o 

Next, introduce x = S/S and set V(x, t) = V/S. This reduces the dimension of 
the PDE by one 

' 9t + ±a 2 (x-<t>) 2 di\v = 



The resulting PDE is closely related to the one found by Rogers and Shi | RS95 1 ■ 
In fact, they can be transformed into each other by a variable change y = x — 4>. 
But our derivation of the PDE seems to be more natural: working in a basis 
of tradables guides us in the right direction. At this point, we will make the 
specific choice of an exponential weight function 

e -i{T-t) 
w(t) = y 

Remember that if the interest rate is constant and equal to r, then the choice 
7 = r leads to an equally weighted average in terms of the dollar price of the 
stock. We now find 

1 - e -^ T -*) 

m = —^r- 

A change of variables is now in place 

This transforms the PDE to 

-d s + ((k + l)z - l -z 2 ^j d z + l -z 2 d 2 )j V = 
A Laplace-transform with respect to s yields 

-A + Uk + l)z - l -z 2 \ 3 Z + \z 2 d 2 \ u = -f(z) 



2 
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where u(z, A) is the transformed function, and f(z) — V(z, T) denotes the payoff 
at maturity. Next, let us define 



u = e? z z K 1 w 



Then w satisfies 



^ - 1 + + ) w = -2e-"z"-V(«) (9) 



where we introduced 



2; 

The homogeneous part of Eq. ^ is Whittaker's equation, and its solutions are 
the Whittaker functions M ti+1 , ll (z) and W K+ i 4l (z) (see appendix). To solve 
Eq. H we will make use of a Green's function approach. A Green's function with 
proper behaviour at the boundaries is given by 

G(x,y) = ^(M K+ i, l t(x)W K+ i, l t(y)l x <y + W K+ i^{x)M K+x ^(y)l x> ^\ 
where Q is the Wronskian of the two solutions 

r(i + 2n) 



W K+ht ,(z)d z M K+ i^{z) - M K+hll {z)d z W K+ i^{z) 



In terms of this, the solution can be written as 

/>oo 

u(z,X) = 2e^z- R - 1 / G(x,z)e-i x x K - 1 f(x)dx 
Jo 

The ASP payoff defined in Eq. @ corresponds to the choice 



f(*)=[^r-*) =k[--l 



a 2 Tz ) \z ) a 2 kT 

Inserting this in the integral, we find that the solution falls apart in two ranges, 
z < a and z > a (or, equivalently, e~ 1T kS < AP and e~ 1T kS > AP). For 
details of this calculation we refer to the appendix. In the former case, we find 

, x 2fce i ^z- K - 1 a K T(-i - k + 

u(z,x) = - ^ y-w K -i,M M *+iA*) 

2k[z{{n - ±) 2 - m 2 ) + a(z - ((« + i) 2 - ,, 2 ))) 

+ z((«-i) 2 -/i 2 )((^+i) 2 -M 2 ) 

The second term is exactly the Laplace transform of x — k. In the latter case 
we find 

, „ 2fce £ ^V K - 1 a'T(-± -K + p) 
(— 2 + K + m)U + « + A*)r(l + 2//) 
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Therefore, the solution can be written as 

v ASP( fc 'T' i ' T ) = 
where I\ and I2 arc defined by inverse Laplace-transforms 



kS(t)h + S(t) - kS(t) e^ T kS < AP 
kS{t)I 2 e~~< T kS > AP 



,p+*oo r( _ 1 _ K + fl )W K _ 1 ^ a )M K+1 ^(z)e > 



e 2 z 



T(l + 2/i) 

r(-| - k + //)M )I _i, M (a)W tt+ i, M (z)e ; 
(-5 + « + /i)(§ + k + A*)r(l + 2/i) 



-dA 



-dA 



' p—zoo V 2 >'^2 

where p is an arbitrary constant chosen so that the contour of integration lies to 
the right of all singularities in the integrand. These integrals can be evaluated 
numerically [AW95, Sha9S]. Note that the value of an average strike call (ASC) 
follows simply from put-call parity, that is, we use 

(s{T) - kS(T)^J - (kS(T) - S(T)J = S(T) - kS(T) 



and find 



V ASC (k, % t,T) 



kS(t)h e-~< T kS < AP 

kS(t)I 2 - S(t) + kS(t) e-~< T kS > AP 



The expression for the value of the ASP at t = 0, i.e. its unseasoned value, 
simplifies considerably. In this case z — *• 00 and s = <j 2 T = 2/(ak), and we find 



V ASP (k n ,0,T) = kS(0)- 



e -§a K f +iCO r(-§-K + /i)M„_i iM (a)exp(2£) 



TV! 



p—too 



(-i + k + + k + A*)r(i + 2^) 



dA 



3.6 T-duality and average price options 



In this section we will use the T-duality, found in section 3.3, to relate prices of 
unseasoned average strike and average price options. We will focus on the case 
where the weight function is exponential. To indicate the weight function used 
in the definition of S, we use a subscript 7, i.e. we define 

Sy{t) = ~ f\-^ T - s ^Y s {t)ds 

It is a straightforward calculation to see that the action of the duality Eq. |(] 
(remark: we use a continuum limit of this result) on the tradables S* 7 , S and P 
is given by 



5 7 (T) «-> e^ T S^(T) 
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Applying this to the payoff of an average price call gives 



KP{0) „,„A + 



(S 7 (T) - KP{T))+ ~ e-^ \S^(T) - — ^- <?(T) 

i.e., it transforms it into the payoff of an ASP. Therefore we see that the unsea- 
soned value of an APC can be expressed as 

V APC (K, 1 ,0,T) = e-^V AS p ( e ^rg( ) 0) '-7.0,r) 
Similarly, we obtain the value of an unseasoned average price put (APP) 
V AFF (K, 7 , 0, T) = e^ T V ASC ( e K ,Ts(oy ~ 7 ' °' T ) 



In this way, we actually reproduce the well known results by Geman and Yor [GY93] . 
3.7 On seasoned Asians 

It is a well-known fact that the price of a seasoned average price option can 
be expressed in terms of the price of an unseasoned average price option with 
a different strike. Let us look at the mechanism behind this. We consider an 
exponentially weighted, continuously sampled Asian with a total lifetime of M, 
expiring at time T (so it is initiated at time T — M < 0) and we are interested 
in its price at t = 0. As before, the payoff of such an option can be expressed 
in terms of tradables S, P and 



T—M 

where we explicitly show the longer sample period in the definition by the sub- 
script M. Now if t £ [0, T] we can write 

s ^ = Tj ir + ^f ^-•'f^ ^ (t) + AP(t) do) 

M M Jx—M P{ s ) M 

where A is proportional to the average over the time period up to t = 0. Sub- 
stituting this in the payoff of an APC, we get 

5 7 , M (T) - KP(T)Y = ^ (s 7 (T) - kP(T)Y 
with 

r> MK - A 

K = 

T 

This shows that the value of the seasoned APC that we are considering can be 
expressed in terms of the value of an unseasoned APC with a modified strike as 
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Of course, the same trick also works for average price puts. Note that it is 
possible for the strike K to become negative. In that case, the option becomes 
trivial. One might wonder what happens if we substitute Eq. [l^ in the payoff 
of an ASC. It turns out that in that case, things do not combine in a nice way. 
Indeed, we find 

(km - s,mt)Y - ( Mts{t) - TS f ] - Ap(T) y (id 

Fortunately, we already have an expression for the value of a seasoned ASC. So 
we can use the formula in reverse, to price options with a payoff given by the 
RHS of Eq. [n]. This will turn out to be useful in the next section. 



4 Cash-dividend 

It well known in the literature how to price options on a stock paying a known 
dividend yield. However, in many cases it is more realistic to assume that the 
cash amount of the dividend rather than the yield is known in advance. This 
makes the pricing problem considerably harder. In this section we show that 
the problem is equivalent to the pricing of Asian options. In fact, we show that 



their prices are connected by the put-call symmetries of section 2.3. This allows 
us to use all the techniques for the valuation of Asian options in the context 
of options involving cash-dividend. The setting is as follows. We assume that 
the stock follows geometric Brownian motion between dividend payments, with 
fixed volatility a (taking the bond as numeraire). Dividends are paid at a set of 
discrete times {ti}, i — 1,2,3,..., < t\ <t% < ■ ■ ■ and are expressed in units 
of the bond 5(U)P. Since we will be interested in options with maturity T, we 
will use a bond with this same maturity. We assume that the interest rate is 
fixed and equal to r, so the bond has dollar value e~ r ( T ~ t \ By Si(t) we mean 
the price of the stock between U and ii+i, in other words 

S(t) = Si{t), for U < t < t i+ i 

So at t = ti a portfolio consisting of 1 stock becomes 

S Q ^6(t 1 )P + S 1 

In order to avoid arbitrage, we assume that the left-hand side equals the right- 
hand side at t\. This can be used to extend the definition of So to all t < t% as 
follows 



S (t), h<t<t 2 



Note that So, by construction, does not make a jump at t\. In fact, So cor- 
responds to the value of the self-financing portfolio that one gets by directly 
reinvesting the cash-dividend payment into the stock again. So So is a tradable 
object. We can repeat the process for the dividend payment at t 2 

S 1 ^S(h)P + S 2 
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Again, S 2 can be expressed in terms of So, extending the definition of So to all 

t < h 

s(t 2 )p(t 2 y 



S(t) = S 2{t ) = (l-^^)S l{ t) = 
1 6( tl )P( tl ) 5(t 2 )P(t 2 ) \ 

By repeating this process, we find that the value of a portfolio V which we 
get by starting with one stock at t = and holding it, together with all its 
cumulative dividends up to time t (note that this portfolio is also a tradable 
object, while the stock by itself is not) is given by 

V(t) = S(t) + J2 S(u)P(t) = ( 1 - E -%7rrM So(t) + £ S(u)P(t) 

where Sq just follows a lognormal price process. Now if we consider a European 
option of the stock with maturity T, we can define the cumulative dividends up 
to maturity as follows 

C(t) = S(ti)P(t) 

Using this we can write 

V{t) = S (t) + C(t)-P(t), t<T 

where 

t<U<T U <t So( - U > 

In terms of these new tradables, we can write 

S(T) = S (T) - P(T) 

Now the connection with Asian option becomes clear. They can be transformed 
into each other by the exchange of So and P, i.e. by using put-call symmetry. 
In fact, we can introduce tradable objects, similar to the Y s (i), as follows 



X s (t) 



In terms of these, we can write 



P(t) t < s 

g&So® t> s 



P(t)=^S( ti )X ti (t) 

t f <T 

Again, for t e [0, T], we have 

P(t) = X T {t) 

So(t) - §^Mt) 
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Therefore we see that the payoff of a plain vanilla option on a stock paying cash 
dividends takes the general form 

with certain weights S. Taking So(t) as numeraire, we see that the tradables 
X s (t) satisfy 

dX s {i) = l t<s aX s (t)dW(t) + --- 

and we can use the integral approach described in section [T^ to price the option. 
Alternatively we can use a PDE approach, i.e. we generalize the definition of P 
as follows, cf. the steps in 3.4, 

P{t) = { 5{s)X s {t)ds = <p(t)P(t) + A(t)S (t) 



A PDE which describes the price process of an option depending on5o, P and 
P (this class includes plain vanilla options on stocks paying cash-dividends) can 
now easily be derived (just take the stock as numeraire) 

dt + ^cr 2 P 2 (dp + cpd P ) 2 ^V = 

It is instructive to consider a call option on a stock which pays a continuous 
stream of cash dividends, with exponential weights. One can use this as an 
approximation to the value of an option where the underlying stock pays a long 
stream of discrete cash dividends. So let us define 

P 1 (t)^^j\-^ T ^X s (t)d S 

The natural choice is 7 = — r, which corresponds to a constant dividend stream 
in terms of dollars. The payoff of a plain vanilla call becomes 

S(T) - KP(T) ) ( S - 8TP- r (T) - KP(T) 



where 5 parametrizes the dividend stream. We will from now on omit the 
subscript from the Sq. By exchanging S and P, exploiting put-call symmetry, 
we get 

P(T) - 5TS^ r {T) - KS(T)^J 

which is the payoff of some Asian option. In fact, by using T-duality, the payoff 
can be related to one that corresponds to a seasoned average strike call (see 
Eq-0) 

{W) S{T) " STerT ^ T) - W) KP{T) ) + 

and we can use the analytical results that we derived for this type of option to 
write down the price of this instrument. 
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5 Conclusion and outlook 



In this article we have shown the power of symmetries to derive prices of com- 
plex exotic options. We focused on arithmetic average options in a Black-Scholes 
setting. By choosing an appropriate basis of tradables, i.e. self-financing port- 
folios, it becomes a straightforward matter to write down the governing PDE 
for the option price. We then proceed to derive the Laplace-transformed price 
of a Euro pean a verage strike option. This result extends the result of Geman 
and Yor [ GY93] for the average price option. Next we show the power of the 



underlying symmetry, by showing the equivalence of the unseasoned arithmetic 
average strike and price options after a suitable transformation of parameters. 
Seasoned options can be treated in a similar way. Finally we exploit the symme- 
try in the problem to show that vanilla options on stocks paying cash dividends 
are equivalent, after suitable transformations, to arithmetic Asian options, thus 
providing a method to price these type of options. 

Let us remark that the present discussion carries over without too much changes 
to the case of basket options and swaptions. We will discuss this in future work. 
Also we did not discuss the case of an arithmetic Asian option with early-exercise 
features. This is however simple to implement and we will come back to this in 
a future work. 



A Whittaker functions 



In this appendix we enumerate some useful properti es of Whittaker functions. 
More information can be found in e.g. [ AS64 , [PBM86 1. The Whittaker functions 
M Ktll (z) and W K ,n(z) are solutions to Whittaker's PDE 



2 _ 1 K (j-M 2 ) 

z 4 z z 2 



f = 



These functions are defined as 

M K ,p(z) = e-5*a" + *iPi(| + n - K, 1 + 2fi,z) 
where the confluent hypergeometric function is given by 

r(6) ^ T(a + n)z n 



L Fi(a,6, z) 



T[a) ^ Li/. -f 



and 



W K Jz) = e~* V+**(| +n-K,l + 2n,z) 
where "L is the Tricomi function 

1 



tf(a,M) 



r(a) 



fc-a-l 



dt 



o 
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There are the following interesting relations 



r'-'M„(i) = H-'-iM,,)-;] 
To evaluate the price of the ASP we made use of the following definite integrals 

J e-^z K - x W K+x ^z)dz = -e-i*z K W K ^z) 

J e-^z K - 2 W K+h ,(z)dz = e-i z z K -\W K -^(z) - W^(z)) 

J e-^z K - x M K+x ^{z)dz = r( 3 + 2 ^^ r(i + k + n)M Ktll (z) 

Je-^z^M K+lA z)dz = ^0^- )X 

x (T(| + k + n)M Ktli {z) + r(-i + k + n)M K _i ifl (z)) 
To calculate its unseasoned value, we recall the asymptotic behaviour of W Kili (z) 
W Kill (z) ~ e~* z z K , z^oo 
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